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OVERSET GRIDS AND HOLE-CUTTING 
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Offset hole 

Create appropriate offset from wall so that interpolation occurs away 
from high gradient regions near wall 


Field equations not solved at blanked points 


OVERVIEW 
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- Brief review of overset grids and hole-cutting 

- Automatic hole cutter closure 

- Adaptive X-rays 

- Automatic hole boundary adjustment 

- Comparisons, test cases 

- Summary and conclusions 


VB DOMAIN CONNECTIVITY / OVERSET ASSEMBLY® 

si mes Research Center 

Hole Points Identification Methods (partial list) 

- Cartesian hole map 

- Line of sight 

- X-rays 

- Implicit hole-cutting 
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Stencil Search Methods (partial list) 

- Cartesian inverse maps 

- Multi-level spatial partitions (ADT, octree, ...) 


Software (partial list) 
- PEGASUS5 

- OVERFLOW/DCF 
- SUGGAR++ 

- PUNDIT 

- OVERTURE 
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ORIGINAL X-RAYS METHOD (OVERFLOW/DCF) 


Meakin, AIAA Paper 2001-2537 


z 

A 



Inside 


Outside 


Rays from image plane 


Pierce points 


□ Point from 
another grid 


Image plane 
spacing As 





APPLICATIONS USING X-RAYS APPROACH ^ 

and Liquid Rocket QM^yctomcj Rotor Blades in Hover 

-77 million “pts i ildT*' ^ Baseline 

JtV 9ridS ^ r " 35 million P ts 

'*£ 1 fHJJ| & " 58 9 r ' ds 

is Z Mm IP ^ Adapted 

“S ° B \ r - 660 million pts 

O | Hi ' T -14000 grids 


Launch vehicle stage- 
separation 
224 volume grids 
267 million grid points 
1 42 X-rays ji 


Launch Site Analysis 

129 volume grids 
92 million grid points 
32 X-rays 
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PROS AND CONS OF ORIGINAL X-RAYS APPROACH 
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PROS 

- Fast hole cutting => excellent for relative motion problems 

- Low memory requirements (2-D map + pierce points) 

CONS - manual inputs needed at the start 

- Identify components 

- Close component open boundaries (tedious) 

- Specify grid subsets to be cut by each X-ray 

(tedious and error prone) 

- Specify uniform X-ray image plane spacing 

(very large X-ray files in tight-gap cases) 

- Specify constant offset for each minimum hole 

(variable offset preferred for better interpolation) 
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COMPONENT DEFINITION 




- Typically a geometric part of a configuration 

- Defined by a collection of surfaces that may or may not 
form a closed volume 

- Each grid point in a near-body volume grid is associated 
with a component 

- A component cutter can blank grid points associated with 
other components 


OBJECTIVES 
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Explore methods to: 

- reduce required user expertise, input effort, time 

- automatically generate appropriate variable offset 
hole boundaries 

Maintain algorithmic efficiency of original X-ray scheme 
for relative motion problems 

Develop code in library form for easy interface with other 
software modules 
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COMPONENT SURFACE DEFINITION 
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Defined via solid wall boundary conditions on surface 
subsets (additional tag to indicate component ID) 

or 

Defined via surface grid cell tags from CAD component 
information (future implementation) 
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AUTOMATIC HOLE-CUTTER CLOSURE 
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WING ROOT CLOSURE 

flifc 

'M 




Bounding box set over open boundary 
(used to clip cells from other components) 

Obtained by sub-dividing longest 
direction of open boundary bounding box 



Clipped surface cells from 
fuselage used to provide 
closure (no need to be 
point matched) 
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COMPONENT OPEN BOUNDARY CLOSURE 




Typically needed at 
junctions between 
intersecting components 
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ADAPTIVE X-RAYS 
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PROBLEMS WITH UNIFORM X-RAY MAP IN 
ORIGINAL SCHEME 




(Arsenio Dimanlig) 



Blade length - 1 95 in. < Q 1 
Chord ~10 in. 


Extremely fine X-ray maps are sometimes needed in tight gaps 
=> very large X-ray file size due to uniform spacing requirement 

Blade and flap - 10 X-rays 

- X-ray file size = 185MB (about 1/4 size of grid file) 
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ADAPTIVE SECONDARY X-RAY MAP 

Proximity in Z Direction 




Abort 


Refined cells (proximity 




Given component primary X-ray map, locate cells in image 
plane for surface grid points P from all other components 


Mark cell for refinement if z-coordinate of P is within s of z- 
coordinates of any pierce points from rays bounding the cell 
(s - As of primary X-ray map) 


4 


Ames Research Center 

Pierce points 
on body \ 



ADAPTIVE SECONDARY X-RAY MAP 

Proximity in X and Y Directions 


Rays from 
image plane 
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Uniformprimary 
image plane 


Mark up cells where 
ray from at least one 4 
vertex has no pierce # 
points 



r fife? 

Each refined cell | ! - 
contains a secondary!^ 

X-ray map 
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ADJUSTED HOLE BOUNDARIES FOR ROTOR 


FLAP TIGHT GAP PROBLEM (SIDE VIEW) 


X-rays file size 
Original = 186 MB 
Adaptive = 63 MB 
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ADJUSTED HOLE BOUNDARIES FOR ABORT 
MOTOR AND MAIN BODY PROXIMITY PROBLEM 




X-rays file size 
Original = 15.3 MB 
Adaptive = 2.9 MB 
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AUTOMATIC HOLE BOUNDARY ADJUSTMENT 
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ACCEPTABLE HOLE BOUNDARIES 




Ideal scenario 

Interpolation occurs between 
cells of comparable attributes 
(cell volume, aspect ratio, 
orientation) 

X-rays approach 
Interpolation occurs away 
from solid wall boundaries 
(high gradient regions) 

Many acceptable solutions 


Aero loads sensitivity? 
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TYPES OF HOLE BOUNDARIES 
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Disjoint Near-body / Near-body 


Near-body / Off-body 


Intersecting Near-body / 
Near-body 
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STEP 1: HOLE BOUNDARY ESTIMATE USING 
WALL-DISTANCE FUNCTION 


Objective: get hole boundary away from minimum hole 
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Standard wall distance function Dw computed in flow solvers with 
popular turbulence models (SA, k-oo, SST) 

= distance from a grid point to closest wall from any component 



Standard wall-distance function 



Hole boundary first estimate 
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HOLE BOUNDARY ADJUSTMENT STRATEGY 



Step 1 

Use a wall-distance function to get first estimate of hole 
boundary location away from minimum hole 

Step 2 

Iterate on hole boundary using orphan points removal as 
objective 



Minimum hole - need to 
expand hole but standard 
wall-distance function is 
zero everywhere on surface 


Need a wall-distance 
function that is sensitive to 
distance to all components 
except self 



DUAL WALL-DISTANCE FUNCTION 
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Dn - Closest distance from point to wall of associated 
component 

- Defined for near-body grids only 

- Distance along ray from wall 


Df - Closest distance from point to wall of any 
unassociated component 

- Defined for near and off-body grids 

- Need to return the ID of the closest 
component (Cn) 

- An approximation may be sufficient 
(pixel map analogy algorithm) 
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STEP 2: HOLE BOUNDARY ITERATION TO 
REMOVE ORPHAN POINTS 


Objective: reduce number of orphan points to same 
number as that in minimum hole 
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- Store interpolation stencils for fringe points visited 
to minimize search at each iteration 


Orphan point in interior 

Convert neighboring blanked points to field 


Orphan point on grid boundary 
Convert blanked points to field 
points near donor grid cell 

Field point conversion not 
allowed at minimum hole points 
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HOLE BOUNDARY OFFSET ESTIMATE USING 
DUAL WALL-DISTANCE FUNCTION 
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Blank point if it is less than about half way from component surface to 

- component grid outer boundary 

- collar grid outer boundary 
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_ CPU TIMES (PERCENT OF TOTAL) FOR VARIOUS 
““'"“§ TEPS OF DOMAIN CONNECTIVITY PROCESS 
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(Not fully balanced OpenMP) 


% of 80 
total 70 
CPU 
time 60 

50 

40 

30 

20 

10 

0 



4.4m 9.9m 16m 28.6m 77.7m 



■X-ray Creation 

■Dual Wall Distance 
Function 
Stencil Search 

■Hole Boundary 
Iteration 
“Other 

(i/o, hole-cutting, etc.) 


(# Grid points, m=10 6 ) 
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LAUNCH VEHICLE (LV) 



Core stage , RSRB’s, MPCV, abort motors, forward and aft attach 

hardware 



(7 million pts) C3P execution time 
25 components = 15.3 min (8 procs) 
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TEST CASES 



- Current software 

- Chimera Components Connectivity Library (C3LIB) 

- Chimera Components Connectivity Program (C3P) 

- Check variable offset hole boundary quality 

- fringe points away from walls 

- Comparison of current scheme (C3P) with 
original X-rays (OVERFLOW/DCF) 

- aerodynamic loads 

- CPU time 

- user time 
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LOCKHEED MARTIN N+2 SUPERSONIC MODEL 
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Fuselage 

Wing 

Blade mount 

Lower wing diverter/nacelle 
Upper wing pylon/nacelle, vertical-tail 

70 near-body grids 
(15 million pts) 

1 off-body grid 
(37 million pts) 

1 component 


C3P execution time 
= 5.7 min (8 procs) 






Lift coef. 

Drag coef. 

Pitch moment coef. 

Original X-rays 

0.4560 

0.0358 

-0.02898 

Improved X-rays 

0.4625 (1.4% diff) 

0.0360 (0.6% diff) 

-0.02967 (2.3% diff) 
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DRAG PREDICTION WORKSHOP (DPW4) 
17 grids, 16.8 million grid points 
Comparison between OVERFLOW/DCF and C3P 





Lift coef. 

Drag coef. 

Pitch moment coef. 

Original X-rays 

0.4852 

0.02741 

-0.02898 

Improved X-rays 

0.4862 (0.2% diff) 

0.02736 (0.2% diff) 

-0.02967 (2.3% diff) 
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DRAG PREDICTION WORKSHOP (DPW4) 
17 grids, 16.8 million grid points 
Various grid slices from C3P 
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CPU AND USER TIME COMPARISON 



CPU TIME OVERFLOW/DCF (original) - 8 MPI processes 

C3P (improve) - 8 OpenMP threads 


Test Case 

# Grid pts 
(xIO 6 ) 

OVERFLOW 

(Standard Distance Fun.) 

C3P 

(Dual Distance Fun.) 

DPW4 

16.8 

1.1 min. 

1.75 min. (+54%) 

D8 

77.7 

14 min. 

19.9 min. (+42%) 


New code is ~ 50% more expensive compared to original 
but it is also doing more work 


USER INPUT PREPARATION TIME 


Test 

Case 

# Grid pts 
(xIO 6 ) 

# 

Grids 

# 

X-rays 

OVERFLOW 

C3P 

Speed up 

DPW4 

16.8 

17 

3 

30-120 min. 

10 - 20 min. 

X 

CD 

X 

CO 

D8 

77.7 

30 

4 

0.5 - 1.5 days 

20 - 30 min. 

36x - 72x 

LV 

65 

67 

25 

1-3 days 

45 - 90 min. 

32x - 48x 
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UNSTEADY 2-D HIGH-LIFT SYSTEM 
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SUMMARY AND CONCLUSIONS 

Minimum Hole Automation 

- Hole-cutter open boundary closure 

- Determination of grid points to be considered for blanking 

- Adaptive X-rays to handle components in close proximity economically 

Adjusted Hole Automation 

- Dual wall distance function to get first estimate 

- Further adjustments using orphan points removal 

- Spatially variable offset for relative motion problems 

Comparison with Original X-rays 

- Preliminary tests show aero loads are comparable 

- CPU time is about 50% more expensive (could be further reduced) 

- Human effort, time, and expertise reduced significantly 
(factor of - 3 - 70) 




Spatially variable offset during relative motion simulation 
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BACKUPS 




r i unc wums. 


- Implement multiple entry points into C3P 

- Improve hole boundary iterations to get orphan points 
back to same number as that for minimum hole 

- Improve optimization and load balancing further 

- Dual wall-distance function 

- Stencil search 
-MPI 

- Systematic study on sensitivity of aerodynamic loads 
to hole boundary locations 
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MINIMUM HOLE CUT USING ADAPTIVE X-RAYS 


z 

t \ 



inside 


outside 


Create secondary X-ray 
map in each refined cell 


Given test point (xp,yp,zp) 

1. Look up cell in primary X-ray using xp, yp 

2. If not refined cell, do inside/outside test on primary map 

3. If refined cell, do inside/outside test in cell secondary X-ray map 
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OVERSET GRIDS AND HOLE-CUTTING 




Minimum hole 
Identification of 
grid points inside 
solid boundaries 

Offset hole 
Create appropriate 
offset from wall so 
that interpolation 
occurs away from 
high gradient 
regions near wall 

Field equations 
not solved at 
blanked points 
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_ APPROXIMATE COMPUTATION OF DUAL WALL- 
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DISTANCE (pixel map /compound eye analogy) 
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Level-1 3-D Cartesian 
map over component 

Level-2 3-D Cartesian 
map over cells with 
walls 

Cw = closest cell 
containing a wall 
located by radial 
search in index 
space 

Df = distance from pt 
to centroid of Cw 
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HOLE BOUNDARY OFFSET ESTIMATE USING 
DUAL WALL-DISTANCE FUNCTION 



Given point from another component, blank point if 
Df < Rb x Db Df < Rb x Dc 

Rb = 0.5 (can be modified by user if needed) 
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